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In this comment, we would like to point out possible 
numerical errors in a recent Letter by Deych et al. jjj . 

In 0, the authors pointed out an interesting ob- 
servation that in some ID random periodic-on-average 
systems, states from pass bands and band edges of the 
underlying band structure follow single parameter scal- 
ing with universal behavior, whereas states from the 
interior of the band gaps do not have universal behav- 
ior and require two parameters to describe their scal- 
ing properties. When the degree of disorder exceeds 
a certain threshold the single parameter scaling is re- 
stored for an entire band gap. Upon careful inspection, 
we find that while these conclusions are qualitatively 
valid, there may be some numerical errors in the results 
of 0. 
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FIG. 1. Band structure 



The model and all the parameters used in the fol- 
lowing computation are the same as in |Q. The model 
in Q is essentially the same as the second model in 
First we compute the band structure using Eq. (6) 
of Q the result is presented in Fig. 1. It shows that 
the first band gap occurs between k = 1.456 and 1.543 
and there is a very narrow second gap located between 
k = 2.992 and 3.005. The third band gap is at around 
k = 4.5. In Q, however, the results indicates a band 
gap at around k = 2.86 which is absent in our compu- 
tation. To validate our program, we used the parame- 
ters from |^| in the computation and obtained the same 
band structure results as in j|. Moreover, by adjusting 
the parameters we were also able to reproduce almost 
exactly the same result for the disorder dependence of 
the Lyapunov exponent as in Fig. 1 of ||. 

Following 0, we then compute the disorder depen- 
dence of the Lyapunov exponent on the standard devia- 
tion in the randomness of the layers' thickness denoted 
by a, and the frequency dependence of the Lyapunov 
exponent and its variance for frequencies covering the 
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first three band gaps. The results are shown in Fig. 2. 
It is clear that the results for frequencies covering the 
first band gaps are qualitatively in agreement with the 
results in Fig. 1 and Fig. 2 of M, except that the band 
gap in [0 is implied at around k = 2.86. In partic- 
ular, as the randomness increases, all the Lyapunov 
exponent curves merge. For frequencies in the vicinity 
of the second band gap, however, the exponent curves 
come across into each other at a certain level of disor- 
der, but split again as the amount of disorder increases 
further. This is shown in Fig. 2(a2). In this case, for 
the same level of the disorder the double peaked struc- 
ture of the variance observed in the first band gap dis- 
appears (Fig. I(b2)). When the disorder decreases, to 
a = 0.00577 say, the double peaks appear again. For 
frequencies near the third band gap, as indicated by 
Fig. 2(a3) and (b3) the results appear to follow the 
features at frequencies around the first band gap. For 
comparison, we also plot in Fig. 3 the dependence of 
the Lyapunov exponent on the standard deviation a 
of the layer randomness at the frequencies indicated in 
Fig. 2 of At these frequencies, all curves seem to 
merge. 
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FIG. 2. Left panel: The dependence of the Lyapunov 
exponent on for different k in the around the first, second 
and third band gaps respectively. Right panel: Frequency 
dependence of the Lyapunov exponent and its variance for 
the frequencies covering the first, second and the third band 
gaps. For the right panel, a is taken as 0.0577. 
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FIG. 3. The dependence of the Lyapunov exponent on 
the standard deviation a of the random layers' thickness 
for different k. 
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